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We propose a laser cooling concept for translational motion of neutral molecules. The cooling 
works by repetition of a cooling cycle composed of three sequential steps: velocity selection of a 
narrow momentum width, deceleration of velocity selected molecules and irreversible accumulation. 
Those molecules which have been cooled undergo only one optical spontaneous emission for irre- 

rf\ ■ versible accumulation. This circumvents the problem of broken closed pumping cycle required for 

laser cooling by optical spontaneous emissions. We consider thermalized molecules at temperature 
of IK as precursor and describe a rotational pumping scheme to prepare the molecules in an ap- 
propriate initial states for translational cooling. Using realistic cooling parameters, we estimate the 
(-H ' cooling time and show that efficient cooling to a momentum width of ~ hk can be achieved within 

a practical timescale. The physical feasibility of the cooling concept is showned by developing a 
cooling model. This enables analytical description of the transient populations and entropies. Sim- 
ulation shows that after the cooling process, a large momentum width is reduced to a narrow final 
width, determined by velocity selection. Translational cooling occurs through the reduction of the 

■^1- ■ center of mass entropy, mainly during coherent laser interactions. The cooling process removes the 

£> ' translational entropy irreversibly to the radiation and internal entropy sinks. 
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I. INTRODUCTION 



The existing laser cooling schemes for atoms have been successfully based on the concept of repeated optical 
excitation followed by spontaneous emission [1]. In molecules, there are additional motional degree of freedoms, 
pH . namely the vibrational and rotational motions in addition to the center of mass translational motion. This creates 
Q_|| two characteristics in molecules which are absent in atoms. First, the optical spontaneous emission (fluorescence) will 
inevitably repopulate more than one state. This makes existing cooling schemes for atoms unapplicable to molecules. 
Second, the existence of additional quantum states associated with the vibrational and rotational motions provide 
additional entropy sinks other than the radiation reservoir for translational cooling. 

A scheme for sequential translational- vibrational-rotational cooling [2] has been proposed but not demonstrated yet. 
Existing schemes for sole vibrational cooling [3] use optimally controlled ultrashort pulses. A translational cooling 
scheme which uses cavity for incoherent processes of dissipation has been proposed [4] . Laser cooling of translational 
motion of molecules is more challenging due to the requirement of large momentum kick from optical photons to 
slow a molecule and the difficulty of maintaining an efficient closed optical pumping cycle, since fluorescence from 
the excited electronic state will populate other molecular levels and leads to inefficient returning to the original 
ground state. Closed cycles in far-infrared and microwave regimes are available but may not be favorable due to their 
low spontaneous emission rates. Thus, a closed optical pumping cycle which involves repeated optical spontaneous 
emissions cannot be simply established in molecules for continuous removal of the translational kinetic energy and 
entropy. 

In this work, we circumvent the problem by first establishing efficient closed optical pumping cycles based on coherent 
processes of Raman velocity selection and STIRAP (Stimulated Raman Adiabatic Passage) [6] to continuously reduce 
the translational kinetic energy of a narrow velocity class of molecules to near zero. Then, this cooled class of 
molecules is pumped to irreversibly populate other rovibrational states through a single optical spontaneous emission. 
These molecules need no further recycling and are decoupled from lasers. This is the main ingredient for successful 
translational cooling of molecules. Our cooling scheme do not introduce permanent vibrational heating, since the 
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molecules subsequently undergo infrared decay mainly to the low vibrational levels in about Is. As in existing 
vibrational cooling schemes [3], the decay inevitably leads to some rotational heating. If necessary, this would need a 
separately rotational cooling scheme, which will be reported in our future work. Our aim here is for translationally 
cooled molecules, which would be invaluable assets for molecular physics since they can be confined in shallower and 
more handy traps for convenient manipulations, fundamental studies and further internal coolings. 

In Section II, we elaborate our cooling concept in a general way. Then, in Section III we describe how the molecules 
to be translationally cooled, can be appropriately prepared in practice using a rotational pumping scheme and buffer 
gas cooled molecules [7]. The cooling scheme is shown to be feasible in two aspects; in practice (as shown in Section IV 
by estimates of physical quantities) and in principle (as shown in Section V by modelling of the cooling process). We 
show that the velocity selection, normally used for subrecoil cooling [5] [8] can be practically used to cool molecules 
from IK to Doppler temperature. We model the cooling process with analytical density matrix and study the entropy 
change of the system. The translational entropy is found to reduce mostly during the coherent deceleration and velocity 
selection steps. Although an irreversible process from the single spontaneous emission per molecule is necessary for 
cooling and leads to the increase in the radiation entropy, it does not contribute substantially to the translational 
entropy reduction. 

II. COOLING CONCEPT 

In this paper, we propose a laser cooling concept for the translational motion of neutral molecules. We describe the 
cooling process with an analytical model and demonstrate the cooling results by numerical simulations. The cooling 
process is based on the repetition of a cooling cycle, each composed of three different physical steps: 1) velocity 
selection, 2) deceleration and 3) accumulation. We refer to the molecules as a total system (subscript 'tot') with two 
subsystems: the electronic, vibrational and rotational states as the internal degree of freedom (subscript T) and the 
molecular center of mass as the external degree of freedom (subscript 'cm'). The velocity selection and deceleration 
involve 4 internal levels, 3 ground states \g0, ±) corresponding to the states with magnetic quantum number M = 0, ±1 
and an excited state |e)(M — 0)(Fig. 2a and 2b). The states |e), \g0) and \g+) are used for velocity selection while 
the deceleration are based on the states |e), \g— ) and |<?+). For convenience, we may refer to \g0) and |<?±) as initial 
state and Raman states respectively. 

The translational cooling starts with the first step, which is the velocity selection, (see Fig. 2a and 3a). A fraction 
of molecules with narrow momentum width and mean momentum P is selected from the population in state |<?0) and 
transfered to an empty state \g+) using a 7r-pulse, without populating the excited state. This step is essentially the 
same as described in Ref. [8] except that here we use orthogonal lin — o~ + lasers instead of counter propagating a + — 
o~ lasers, since the dipole moment between states |g0) and |e) is aligned along the a + polarized laser axis. 

In the second step (see Fig. 2b and 3b), the selected molecules are decelerated using a coherent process. In principle, 
non-optical deceleration technique can be applied [10]. However, for practical reasons, we require efficient and fast 
optical momentum transfer by repeated population inversion between the Raman states \g+) and \g—). This can be 
accomplished by various population transfer techniques using optimally engineered laser pulses [11]. Here, we adopt 
the STIRAP process which has been shown to be highly efficient in molecules [6] [12] and has been used for coherent 
momentum transfer in atom optics and atom interferometry [13]. In this step, the internal populations are repeatedly 
inverted between \g+) and \g—) ■ Each inversion is a substep of step 2, which is accomplished by reversing the directions 
and the counter-intuitive sequence of a pair of counter propagating o~+ and a— polarized pulses. Each pair of the 
STIRAP pulse provides about 2Kk of photon momentum kick to decelerate the molecules. Thus, the population with 
a mean momentum P will require a number of P/2hk STIRAP inversions to reduce the mean momentum P towards 
zero. In the end of this step, we have a fraction of molecules with the translational energy close to zero . 

We want to repeat steps 1 and 2 to accumulate the populations with narrow momentum width around zero mo- 
mentum. Since steps 1 and 2 are coherent processes, subsequent repetition of these steps would not lead to net 
cooling because the next STIRAP process would accelerate the molecules out of zero mean momentum. The Raman 
states must be emptied before the next deceleration cycle starts. We achieve this with an irreversible process in the 
third step(see Fig. 2c and 3c). The cooled and decelerated molecules are rapidly transferred from the Raman states 
to a decaying state \d) which may be an excited electronic state or a rovibrational state in the ground electronic 
state. Due to the requirement of fast translational cooling of unconfined molecules, we use an excited electronic state 
which provides fast optical spontaneous decay time. Subsequent spontaneous emissions populate the many vibronic 
ground as the accumulation states {|acc)}. We emphasize that there is no vibrational heating as time passes, since 
the molecules in states {|acc)} are decoupled from the laser beams and will undergo tens of infrared spontaneous 
decays to the lowest vibrational levels in about Is. This introduces additional momentum spread in the order of an 



optical photon momentum, much smaller than the overall reduction of momentum width. It is the irreversibility of 
the spontaneous emission which enables repeated accumulation of the molecules with narrow momentum slices and 
zero mean momentum from each cooling cycle as steps 1, 2 and 3 are repeated. 

Each cooling cycle which is composed of a sequence of the three steps described above, begins with the velocity 
selection of a new slice of narrow momentum width. The whole cooling process involves the repetition of the cooling 
cycle. We point out that the equivalent temperature of the molecules after the process corresponds to the velocity 
selected momentum width, which can be well below the recoil limit. 

The main difficulty to laser cool molecules by repeated spontaneous emissions in the optical regime is due to the 
fluorescence decay from an excited electronic state to many rotational- vibrational states of the ground electronic state. 
After a spontaneous emission, the molecules have a high probability of not returning to the original state. Thus, the 
molecules decay to a state which is off-resonant with the pumping laser before their momenta are reduced by the laser 
photons. In our cooling scheme, the momentum and kinetic energy of each molecule are quickly reduced toward zero 
by deceleration (step 2) before undergoing a single optical spontaneous emission process in the end of a cooling cycle. 
Thus, each cooling cycle completely cools a particular fraction of molecules. This cooling scheme circumvents the 
requirement of repeated spontaneous emissions for each molecule to obtain cooling. It is this unique feature which 
makes translational laser cooling of molecules possible. 

III. PREPARATION OF INITIAL MOLECULES 



Initially, the molecules to be translationally cooled can be prepared in a thermal equilibrium of IK by the versatile 
buffer gas cooling [7], with a single heavily populated internal state \g0) and two empty states |fltt) using a rotational 
pumping scheme described below. 

A. Thermalized molecules 

We start with the preparation of thermalized molecules which can have sufficiently few internal states for convenient 
rotational pumping and sufficiently low translational temperature for practical translational cooling in unconfincd 
space. We emphasize that we are dealing with neutral molecules instead of molecular ions. Neutral molecules 
produced by laser ablationed are highly excited internally and translationally. Buffer gas cooling is a good technique 
to thermalize the molecules, reduce the translational temperature of IK and depopulate many of the highly excited 
vibronic and rotatinal states. Elastic collisions during buffer gas cooling provides a good thermalization mechanism 
to bring the molecules close to thermal equilibrium. At thermal equilibrium of temperature T,, the populations of the 
ground electronic surface |£) are distributed among rovibrational states (n, J) according to the Maxwcll-Boltzmann 
distribution [14] 

P(n, J) ex (2 J + 1) cxp{-hcBJ(J + 1) - hcv(n + -)}/k B Ti (1) 

where v and B are vibrational and rotational constants (in wavenumber) respectively. The maximum population 



corresponds to n max = 0, J max ~ y 2 hcB ~ \- ^ initial temperature of Tj = \K and B w 0.1cm 1 , the molecules are 
mainly in the ground electronic state with the highest populations in n = and J max w 1 or 2. Thus, the molecules 
are reasonably cold internally and this makes it easy to create a single heavily populated state \g0) for subsequent 
translational cooling. Since the excited vibrational states are scarcely populated, molecular collisions rates are small 
and can be neglected [15]. 

B. Rotational Pumping 

We wish to create 2 empty states |<jtt) and a heavily populated state \g0) for more efficient translational cooling. 
This can be done by pumping the populations from the most highly populated rotational sublevels into a single 
sublevel within the ground vibrational and ground electronic state, namely \g0) = |S,n = 0, J = 1, Mj = 0)(see Fig. 
1), Using a + polarized laser, the populations in rovibrational states \n,J,{Mj}} = |0, 1, — 1), |0, 2, {— 2, — 1,0}) are 
excited to states |1, 0), |1, 1, {—1, 0, 1}) respectively. On the other hand, o~ polarized laser (same frequency as er+) 



excites. the populations from states 1 0,1,1), |0,2,{2,1,0}) to states |1,0), |1,1, {1,0,-1}) respectively. A third laser, 
linear polarized with higher frequency excites the population in |0, 0) to |1, 1, 0). Subsequent infrared decays by virtue 
of selection rules, A J = ±1 and AM,/ = 0,±1 effectively transfer the populations from the above states to a single 
state |0, 1, 0)(referred as \g0)). The empty states of |1, 1, ±1) are referred as \g±). The a + and a~ polarized lasers are 
switched on alternatively instead of simultaneously to avoid the creation of a dark state which will prevent complete 
pumping to state |<?0). Taking into account the magnetic levels splitting by the trapping field and small differences in 
the rovibrational transition frequencies, we could use the lasers with a finite bandwidth to resonantly excite all the 
above transitions near the trap center. The pumping process takes about 5tir ~ 0.05s (tir is the IR natural lifetime) 
which is considerably shorter than a typical trapping lifetime (for example CaH lifetime is about 0.6s [7]). 

Even after rotational pumping, there is a small fraction of the molecules remaining in higher vibrational and 
rotational states. These molecules are not in the 4 cooling cycle state (\g0), |</±) and |e)). Therefore, they will not 
be translationally cooled and can be disregarded. Among these 4 states, only a single state \g0) is populated while 
the rest are essentially empty within a rovibrational decay lifetime. Thus, we still have a large fraction of buffer gas 
cooled thermal molecules in a single internal state \g0) as the precursor for translational cooling. When the trapping 
field is switched off, the laser cooling starts. 

IV. ESTIMATION OF COOLING PARAMETERS 

We show that the propose cooling scheme can be realized in practice. Based on the results of our theoretical 
studies, we use realistic physical parameters to briefly estimate the cooling time for each step. It is essential to have 
a short total cooling time for unconfined cooling so that the molecules are cooled before they fly apart beyond the 
laser cooling region. Once translationally cooled, the molecules can be easily manipulated and confined in a simple 
and shallow trap. 

A. Velocity selection 



The velocity selected momentum width can be arbitrarily narrow using long duration Raman 7r— pulses. For 
practical purpose, it is essential to consider an optimum pulse duration for this step. Since the selected population is 
shifted by discrete number of H(k- + k + ) w 2hk in the deceleration step, it is most efficient to have the selected width 
of AP VS rs 2Hk, corresponding to 4 times the recoil temperature, which is T vs = 2 mJ ~ 4^-^ f° r optical wavelength 
of A = 300nm and molecular mass of M = lOOa.m.u. For Blackman pulse shape [9], the required 7r— pulses duration 
is estimated to be n w TEF~ = j^~ ~ 27/is, which is longer than the natural lifetime (r e w 10 7 s _1 ) of state |e) and 
the vibrational period of a molecule [16]. Since the sweep of the Blackman pulses is not a rapid passage, the excited 
state will be appreciably populated unless the optical detuning 5 is large enough, but smaller than the rotational level 
spacing (around 10 10 s~ 1 ) to avoid excitations to untargetcd rotational states. We use 5 — — 50r «5x 10 8 s _1 . The 
corresponding Rabi frequency amplitude is fl = W—Sn/JjiOSV) s=s 9.8 x 10 6 s _1 [17]. 

For a molecular gas of mass lOOa.m.u. and initial temperature Tj = IK, the momentum width is The maximum 
mean momentum is taken to be P max = V^MknTi ~ 1000hk. If each velocity selection step selects a momentum 
width of AP VS = 2hk, the total number of velocity selection steps which is the number of cooling cycle, A max = 
2Pmax/A.P vs w 1000, The total time required for velocity selection in the cooling process is t\ — iV max Ti w 27ms. 

B. Deceleration 



For complete population transfer between states \g+) and \g— ) in each STIRAP inversion, in addition to counterin- 
tuitive sequence and overlapping pulses, the pulses duration T2 must meet the adiabatic criteria [18]. Our simulations 
which take into account spontaneous emissions using master equations give a criteria of T2$l > 50 for 99.95% transfer 
efficiency with zero detuning. With an efficiency of e = 99.95% for each STIRAP inversion, the minimum deceleration 
efficiency corresponds to the maximum number of inversions -/V max = 1000 for molecules with mean P/v max , which is 
e N = € N ^ =60%. 



In one-photon resonant condition to intermediate state |e) = \X, n = 0, J = 0,Mj = 0), the finite spectral width 
Au) of the pulses must be much smaller than rotational level spacing u> ro t to avoid off-resonant coupling to additional 
intermediate states which may lead to imperfect transfer efficiency; for example, Aw rj 1/t2 << uj rot ~ 10 10 s _1 . 
Thus, one can safely use the typical [18] nanosecond STIRAP pulses and have each inversion duration of T2 = 10 _8 s 
corresponding to the spectral width of Alu rj 10 8 s _1 and a moderate peak Rabi frequency fi s=s 5 x 10 9 s _1 . 

The number of STIRAP inversions required to decelerate a velocity selected slice with mean momentum py is 
N s=a Pjv/2fi,fc and the corresponding deceleration time is ijv = Nt2- The total deceleration time in the cooling process 
is estimated as 

N max N max 1 

*2 = 5Z ijv = S Nt2 = 2 iVmax ( 7Vmax + 1 ) 75 - 5ms ^ 2 ) 

JV=1 JV=1 

Despite the large number of STIRAP inversions, the total cooling time in step 2 is still less than the total time for 
velocity selection. 

C. Accumulation 



In practice, the decelerated molecules in \g+) or \g—) can be driven resonantly by c.w. laser to an excited electronic 
state for flourescence decay to many ground rovibrational levels within T3 = 10/r e w 1/is. The total accumulation time 
is £3 ~ -/V max T3 w 1ms. Thus, the total time for the whole cooling process is t coo i = 1 1 + £2 + £3 ~ 33ms. The cooling 
can be done simultaneously on the positive and negative sides of the momentum distribution. If the velocity selection 
and deceleration start from the fastest molecules to the slowest, the drift distance for molecules with momentum P/v 
is estimated as Ln = N2hk{^(N max — N)(N max — N + l)r 2 + (N max — N)(n + t 3 )} which gives a maximum drift 
distance 20cm for N = 475. The maximum drift distance can be reduced by optimizing the order of velocity selection 
on certain velocity class. 

V. MODELLING OF THE COOLING EFFECT 
A. Initial state and density matrices 

From the practical descriptions in Section III, we assume that initially the internal state is pure and uncorrelated 
to the center of mass momentum state. Only the internal state \g0) is occupied. The center of mass subsystem is in 
the thermal Gibbs state [20] and the initial total density operator is written as 

/WO) = /5/(0) p cm (0) = \g0){g0\ ]T W(P)\P) (P\ (3) 

p 

where W(P) = i?e ~ p l a is the statistical weight, Z = J] e~ p l a is the partition function and a — \/2MksT is the 

p 
initial momentum width with effective temperature T. On one hand the gas is coupled to lasers, while on the other 

hand it forms an open system whereby it can exchange entropy with the outside environment. 

The quantum dynamical evolution of the molecules in each cooling step and throughout the cooling process can be 

modeled analytically by using the density matrix elements. This enables convenient evaluation of the internal and 

center of mass subsystems probability distributions and the corresponding entropies, which allow us to work out the 

dynamics and the cooling efficiency. The density matrices of the total system, internal subsystem and center of mass 

subsystem in internal basis {\a}} and center of mass momentum basis {|P)}, are defined respectively as 

/M*) = £X>a6(P,i y ,t)|a,P)<&,P'| (4) 

a.b P,P< 

Pi(t) = Tr cm {p tot (t)} - J2 \a)(b\J2P"b(PPt) (5) 

a, b P 

Pc m (t)=Tn{p tot (t)} = ]T \P)(P'\J2paa(PP',t) (6) 

P,P' a 



where p a b{P, P', t) = (a, P\p to t(t)\b, P') are the total system density matrix elements, a, b £ {e, gO, g±, d, {ace,}} the 
internal states and {P, P'} the center of mass momenta along er± polarized lasers axis (defined as z-axis). 

The time evolution dynamics of the total system in step 1 and step 2 are formally described by the 3-level Bloch 
equations [?]. In step 1 the total system states are \gO,P), \e,P) and \g+,P — hk). A STIRAP inversion from \g+) 
to \g—) in step 2 involves the momentum family of states \g+, P + hk), \e, P) and \g—, P — hk). Due to the reversal 
of the pulses during subsequent inversion from \g— ) to \g+), the family of states becomes \g— , P — hk), |e, P — 2hk) 
and \g+,P — 3hk). With the high populations transfer cfhciency of STIRAP, the center of mass coherences in the 
populations, for example (<?+, P—hk\p tot (t)\g+, P—3hk), are negligibly small. We start with populations in the internal 
state |g0) and the center of mass Gibbs state with no coherences in momentum basis. So, the only non vanishing 
matrix elements are poo(P, Pi 0). The processes of a velocity selection, a STIRAP inversion and a spontaneous emission 
do not create center of mass coherences. Therefore, the populations remain diagonal in momentum basis throughout 
the cooling process: p aa (P, P',t) = p aa (P, P, t)S PtP > . 

On the other hand, transitions between different internal states are accompanied by momentum recoil. Thus, the 
matrix elements between different internal statcs(a ^= b) arc non vanishing only for the center of mass coherences 
P^P 1 : p ab (P, P', t) ^ and p ab (P, P, t) = 0. Accordingly, Eqs. (3) and (4) reduce to 



Pi 



(t)=J2C a (t)\a)(a\ (7) 



Pom(t)=J2f(P,t)\P){P\ (8) 

P 

where the internal and center of mass probability distributions C a (t) and f(P,t) are respectively defined as 

Ca(t) = ^2 P aa(P,P,t) (9) 

P 

f(P,t)=Y,Paa(P,P,t) (10) 

a 

Thus, we only need to evaluate the diagonal density matrix elements p a a{P, P, t) in order to calculate the subsystem 
density matrices and the corresponding entropies. 

VI. POPULATION DYNAMICS 

Step 1 of N-th cooling cycle: To simplify the mathematics, we can model the velocity selected momentum distribution 
by a Gaussian function of a narrow width a vse i(« a). During step 1 of the iV-th cooling cycle (subscript N:l), only 
the populations in the states \gQ) and \g+) change with time and can be written analytically as the following: 

N-l 

Poo(P,P,h) N ..i - w(P) - Y, v j( p ) - v N (P)hi(h) (ii) 

P++ (P, P, h) N:l - h l (t 1 )W(P N )e-( p+ahk »>+- p ^ 2 /^ (12) 

where < ti < T\. t\ is the duration of the population transfer in step 1. The function hi(ti) — sin 2 irt/2Ti 
describes the time evolution of the populations. Vj(P) = W(Pj)e~( p ~ p ^ l a ^^i is the velocity selected distribution in 
the j-th(j = l...iV) cooling cycle with mean momentum Pj, momentum width o~ vse i and mean population W(Pj) = 
\& Pj ' a ■ kN+ is the cr + laser wavevector of the N-th cooling cycle and a = sign(Pj) determines the direction of 
the o~ + laser to provide decelerating momentum kick. 

The first term in Eq. (9) describes the initial populations distribution. The second term corresponds to the depleted 
populations after N—l cooling cycles. The third term describes the transient population depletion in the iV-th cooling 
cycle. In the end of step 1, the populations become 

N 

p 00 (PP,T 1 ) N .. 1 = W(P)-J2V j (P) (13) 

P++ {P, P, T1 ) N:1 - W(P N )e-( p+ahk ^+- p ^ 2 ^< (14) 



Step 2 of N-th cooling cycle: There are many substeps within step 2. Only the populations in states \g+) and \g— ) 
change with time during the n-th substep of step 2 in the iV-th cooling cycle, according to 

Paa(P, P, i 2 )jV:2(n) - h 2 (t 2 )W ' { p N ) e -(P+^ N , ++a K^P N f ', /^„, (15) 

p bb (P,P,h) N:2{n) = (1 - h 2 (t 2 ))W(P N )e-( p+ahk »<++ aA -i -r»)*l°l.« (16) 

p ++ (P,P,0) N:2{1) =p ++ (P,P,r 1 ) N:1 (17) 

where a = — , b = + for odd n(inversion from |.g+) to \g— )) and a = +, 6 = — for even n. The time interval for every 
substep is < t 2 < t 2 , each with the same duration r 2 for complete inversion. h 2 (t 2 ) = sin 2 irt/2T 2 describes the time 

n 

evolution of the populations and A„ = J^ h(kj + + fcj_ ) w 2nhk is the total amount of momentum transfer provided by 

i=i 
n pairs of STIRAP pulses. The precise values of fc/v+, %+ an d fcj- are obtained via the two-photon Raman resonance 
condition. In the numerical simulations they are simply taken to be approximately equal to k without affecting the 
validity of the results. 

The momentum width a vse i is unchanged throughout step 2. The number of deceleration cycles n max required to 
reduce the mean momentum from Pjy to zero can be predicted through aftfc/v,++aA„ — Pn ~ cth~k(l + 2n max ) — Pjq ~ 0, 
which gives 

1 < n < (n max = ^- k -\)- (18) 

In the end of step 2, the populations become 

Paa(P,P,T2) N :2(n^ x ) « W{P N )e~ p2 l^^ (19) 

Step 3 of N-th cooling cycle: After step 2, if the cooled populations from the Raman states |<7±) are transferred 
rapidly to the decaying state \d) with a short 7r-pulse, the process is essentially unitary due to negligible sponta- 
neous emissions from \d). Furthermore, if the laser is orthogonal to the z-direction, no momentum kick is imparted 
in z-direction. This enables the replacement pdd(P, P, 0)jv:3 ~ paa{P,P, T2)jV:2(n max ) as initial population, setting 
p±±{P, P, ta)N:3 = and disregarding the dynamics of the tt — pulse. The population in \d) decay to many vibrational- 
rotational levels in the ground electronic state. Assuming that only a small branching ratio of the population from 
\d) decays back to |<?±), the time evolution of the populations in step 3 is simply due to the spontaneous decays from 
\d) to M number of accumulation states {\accj),j — I...M} with the corresponding decay rates Tj. This is given by 

M 
Pdd(P,P,h)N:3 ^J2 F ^ r3t3W ^ PN ^ P2/al ' 1 ( 20 ) 



IV — J. 

p a cc j {P,P,h)N:^=F j ^^e- p2 /^{Y J W(P 3 ) + W(P N )(l-e- r ^)} (21) 

O ace - t 

where < ^3 < T3 is the time interval for step 3, T3 is the duration for \d) to empty its populations to essentially zero 

M 

(taken as 10/Tj) and Fj is the Franck-Condon factor for level \acCj) with J2 Pj = 1- 

The first term in the curly bracket of Eq. (19) corresponds to the cumulative populations from the previous N — 1 
cooling cycles. It is equal to the depleted populations from \g0) while the second term describes the transient increase 
in populations. The accumulated momentum width o~ acc is slightly larger than a vse i due to the small momentum 
spread (about 0.5ftfc a ) from a single spontaneous emission. The normalization factor 2Wl j n fr^ (19) ensures that 
Trtot{ptot} = 1 and can be estimated by 

O'acc/o'vsel ~ 1 + 0.5hk a / a vse i (22) 

where k a — (E d - E acc )/hc. 

VII. ENTROPIES 

The cooling effect is shown by the variation in the momentum probability distribution and the variations of the 
total, internal and center of mass von Neumann entropies [19]. From Eq. (1) the initial internal, center of mass and 
total entropies are respectively: 



5/(0) = (23) 

S cm (0) = -k B J2w(P)lnW(P) (24) 

p 
Stot(0) = 5 cm (0) (25) 

The expressions for internal and center of mass entropies are obtained from Eqs. (5)-(8) using Sx(t) = 
Trx{px(t) Inpx} where X <G {I, cm}, as follows: 

5/ = -fc B ^C a lnC a (26) 

a 

5 cm = -fc B ^/(P)ln/(P) (27) 

p 

Equations (24) and (25) are valid for all times in all cooling steps and take a simple form because the reduced density 
matrices are diagonal (Eqs. (5) and (6)). The center of mass probability distributions for step 1, 2 and 3 of the AT-th 
cooling cycle are given by: 



f(PM)N:l= J2 Paa(P,P,tl) N :l+J2P^A P > P > T 3)(N-iy.3 (28) 

a=0,+ j 

f{P,t2)N:2(n)= P0o(P,P,Tl)N:l+ ^ Paa{P, P, *2)iV:2(n) + ^2 Pac Cj (P, P, T 3 )(N-1):3 ( 29 ) 

o=+i- 3 

f(P,t3)N:3 = P00(PPT 1 ) N .. 1 +p dd (P,P,t 3 ) N .. 3 + Y / Pacc ] (PP,t 3 ) N .. 3 (30) 

3 

Using J^° e~( p ~ a ) / <7 ««eidP = a vse i^/n, we can show that f(P) is normalized, J f[P)dP = 1. Similarly, the internal 
probability distributions for step 1, 2 and 3 of the iV-th cooling cycle are given by: 

Co,+ (*l)jV:l = / y P0O,++(-P i -P i ^l)jV:l 

p 

Cac Cj (tl)N:l = ^Pacc 3 (P, P,T 3 )(N-1):3 (31) 

P 

Co(t2)N:2(n) = / , POojP, P, Ti)n-.1 
P 

C+,-(*2)jV:2(n) = 2_^ P+ + , {P, P,h)N:2(n) 

P 

Caccj (*2)/V:2(n) = 2_^ Pacc 3 (-P, P, T"3)(jV-l):3 (32) 

P 

C (i3)w:3 = y^/QQo(P,P,ri)jV:l 
P 

Cd(*3)iV : 3 = ^P<ta(-P,-P,*3)j\r:3 

P 

C aC c 3 (t 3 )N: 3 = ^2 PacCj (P, P, *3)jV : 3 (33) 

P 

The iofaZ density matrix in this step involves the basis states \g0), \d) and {|acc.,-)}. It is diagonal since these states 
are not coupled by laser. The evolution is characterized by free spontaneous emission from \d) to {|accj)}.Thus, the 
time evolution of the total density matrix for step 3 can be written in a completely diagonal form 

Ptot(i3)jV : 3 = ^{/900(P,P,Tl)iV:l|0,P)(0,P|+ J2 Paa(P,P,t 3 ) N : 3 \a,P)(a,P\} (34) 

P a=d,{accj} 

with Tr to t{ptot,N(t)} = 1. 



The total entropy changes only during incoherent process in step 3. Its time evolution in the iV-th cooling cycle is 
given by 

Stot(t 3 )N:3 = -k B ^{p 00 {P, P, ri)jV:llnpoo(-Pj P, 7"l)jV:l + ^ Paa{P,P,t3)N:3^Paa{P,P,t 3 )N:3} (35) 

P a—d^{accj} 

Throughout the coherent interactions of step 1 and step 2 in the JV-th cooling cycle, the total entropy takes a constant 
value of 5 , tot (T 3 )( A r_i) : 3. 

Each cooling cycle is completed in step 3, during which radiation entropy Sn(t) is created due to the random 
direction and many frequencies of the spontaneously emitted photons as the cooled molecules decay to the many 
\acc) states. From the knowledge of the angular distribution of the spontaneously emitted photons, we can write the 
radiation density matrix as 



PR 



1 pill pTT 

^-E/ F j M jX (e)\j,x,e, ( (>)(j,\,eA\smededcf ) (36) 



where \j, X,0, <j>) is the state of the photon emitted with frequency u acc . — (Ed~E acc .)/h, polarization A = 0, ±1, and 
in spherical coordinate angles (0, 4>) relative to the dipole axis. M x (0) = |(1 + cos 2 6)5\.±i + | sin 2 0S\,o corresponds 
to the probability of emission at angle relative to the dipole axis and Fj is the Franck-Condon factor. The trace 
over all decay fequcncies and angles gives Trj t x,e,<i>{PR} = 1- The radiation entropy corresponding is easily obtained 

as 

SR(t) = -k B Y] [ R jXe In R jXe d0 (37) 

where R jXd = FjM jX (0) sin 0. 

We note that the existence of many \acc) states during step 3 leads larger radiation entropy due to the frequency 
dispersion, in addition to the angular dispersion of the spontaneously emitted photons. 

VIII. RESULTS AND DISCUSSIONS 

Based on the model discussed above, we simulate the essential features of the cooling process. The typical variations 
of the internal populations distribution for each step in a cooling cycle is shown in Figs. 3. In each cooling cycle, a 
fraction of molecules with a narrow momentum width arc depleted from the state \g0) by velocity selection (Fig. 3a), 
slowed towards zero mean momentum via repeated inversion by STIRAP (Fig. 3b) and dumped irreversibly into the 
accumulation states (Fig. 3c). 

We give the overall results of the cooling process first and discuss the details of the cooling mechanism later. The 
cooling effect is shown by the narrowing of the momentum probability distribution f(P), which evolves from a large 
momentum width (about 30fifc) into a narrow momentum width (about 2hk) with zero mean momentum (Fig. 6). The 
main interest of this paper is the center of mass entropy S cm which is found to decrease after the cooling process (Fig. 
5a), in particular during velocity selection and deceleration since these steps move a fraction of the populations closer 
to zero momentum. This is the sign of cooling in the molecular translational degree of freedom. The coherent processes 
of steps 1 and 2 contribute substantially to the reduction of S cm , from the maximum entropy of Gibbs state to the 
final populations with narrow momentum width. Throughout the cooling, the internal entropy Si increases(Fig. 5a). 
Within each cooling cycle, the laser fields have created correlations between internal and center of mass momentum 
states, thus transferring the entropy from the center of mass to the internal subsystems by reducing the center of 
mass state dispersion while increasing the internal state dispersion. Indefinite repetition of the center of mass cooling 
effect is possible through the irreversible process of step 3. 

We proceed to elaborate on the detailed features of the entropy variations in each cooling cycle with regard to the 
populations of the internal and center of mass states. During velocity selection, the extraction of a narrow slice of 
momentum from \g0) leads to distortion in the momentum distribution. One hk of momentum kick on the velocity 
selected population leads to a narrower overall momentum distribution and to the reduction of S cm (Fig. 4a). At the 
same time, Si increases(Fig. 4b) because an additional internal state \g+) is increasingly populated. 

In the deceleration step, S cm is reduced in each substep (Fig. 4a) by repeated momentum kicks which push the 
selected populations closer to zero mean momentum. Entanglement [21] of the center of mass momentum states in each 



STIRAP process increases the population dispersion of the center of mass states. This explains the transient increase 
in S cm , which is less than the overall reduction of S cm in one STIRAP process. During the first STIRAP process in 
step 2, this transient increase is not observed because it is quenched by the rapid reduction of S cm . The first STIRAP 
process reduces S cm more than the subsequent ones. This quantitative difference is due to the nonlinear dependence of 
the entropy on the momentum distribution. However, Si oscillates above a constant value (Fig. 4b). Each oscillation 
corresponds to a substep, which is an inversion process by STIRAP and can be explained by the creation of adiabatic 
states through the entanglement of states \g+) and \g—}. The entanglement involves the population of an additional 
internal state, thereby increasing Si. After the completion of a STIRAP process, the entanglement disappears and Si 
falls back to the value before the entanglement. Thus, the value of Si before and after step 2 are the same. 

The total entropy S to t is constant throughout steps 1 and 2 but changes only during step 3 (Fig. 4c). This 
shows that unitary evolution can reduce the center of mass entropy despite the invariance of the total entropy. In 
other words, coherent interaction alone is able to cool the center of mass subsystem at the expense of heating the 
internal subsystem. However, the cooling cannot be repeated unless an irreversible process is introduced as part of 
the cooling cycle, which is step 3. In step 3, the slight increase in S cm due to momentum spread from one spontaneous 
emission is smaller than the decrease in S cm in steps 1 and 2 (Fig. 4a). This leads to the overall reduction of S cm 
in each cooling cycle. In fact, it is the coherent process of steps 1 and 2 instead of the spontaneous emission, which 
substantially reduces S cm . This is an essential point, as it indicates that coherent process can lead to substantial change 
in the subsystem entropy. The important role of the single spontaneous emission is to irreversibly accumulate the 
translationally cooled populations. We note that this accumulation role can also be realized if the cooled populations 
in each cooling cycle are coherently transferred to different internal \acc) states. For this, the number of \acc) states 
has to be more than the number of cooling cycles and multi- frequency lasers resonant with all the \acc) states would 
be required. In this case, the whole cooling process is coherent and no spontaneous emission is required. Then, both 
translational and internal cooling cannot be achieved simultaneously. Translational cooling must be accompanied by 
internal heating, or vice versa. Thus, coherent cooling of translational motion is not a conceptual impossibility but 
rather a practical one, at least at the present moment. 

The entropy variations for a single \acc) state and multi \acc) states are shown in Fig. 4b and Fig. 4c (we use 10 
| ace) states with logarithmic Franck-Condon factors and a fixed decay rate). With single |acc) state, Si reduces during 
step 3 of the second cooling cycle because the number of internally populated states reduces from three populated 
states (\g0) , \g±) , \acc)) to two states (\g0), |acc))(Fig. 4b). The center of mass entropy S cm is unaffected by the 
number of |acc) states (Fig. 4a). However, the internal entropy Si increases with more |acc) states since the number 
of populated internal states increases (Fig. 4b). It is interesting to observe that the increase in 5*/ from populating 
10 internal states(in step 3) is still smaller than the increase in Si from velocity selcction(step 1). This is because 
entropy is not simply a measure of the number of occupied states, but also depends on the relative populations among 
the states. 

The internal and center of mass entropies fluctuate throughout the cooling process (Fig. 5a). As the momentum 
width of the populations shrinks towards a narrow velocity selected width in the cooling process, the S cm reduces from 
the maximum entropy state towards zero. The small residual populations in the right and left wings of the narrow 
peak (Fig. 6) is due to inefficient velocity selection with Gaussian momentum profile. In practice, a more efficient 
pulse shape can be tailored to selectively transfer the narrow populations as close to a square momentum profile as 
possible. This would bring the final S cm closer to zero in the end of the cooling process. On overall, Si increases in 
every cooling cycle due to the increase of the populations in the many |acc) states, thus changing from populating an 
internal state (Si = 0) to populating multi internal |acc) states. 

The total entropy Stot before and after the cooling do not seem to change a lot (Fig. 5b). In fact, Stot can increase 
or decrease due to the coupling of the total system to the environment. The total, internal and center of mass entropies 
satisfy the Araki-Lieb inequality [22], |Sj— S cm \ < Stot < Si+ S cm , as shown in Fig. 4c and 5b. Therefore, the 
cooling results are consistent with the general properties of quantum entropy. The non-negligible index of correlation 
Ic = S cm + Si — Stot [20] shows that the internal and center of mass subsystems are correlated by the laser interactions 
throughout the cooling process. However, after cooling the correlation is not maximum since Stot is not zero and 

S C m T 1 Si- 

We note that the increase in the internal entropy is due to the specific choice of the accumulation scheme (step 3), 
due to practical requirement for faster cooling time. The amount of increase in the internal and radiation entropies 
depend essentially on the number of |acc) states. The resulting internal heating is entirely analogous to creation of 
radiation entropy. The many internal states in molecules provide an additional entropy sink, much like the radiation 
Fock states. During the cooling process, the translational entropy is removed irreversibly by spontaneous emissions 
to these entropy sinks. 

If we had coherently accumulated the cooled molecules by transferring each cooled fraction of molecules after each 
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cooling cycle to distinct rovibrational states instead of using spontaneous emission, less radiation entropy will be 
created while the internal entropy increases dramatically. This requires multi-frequency lasers, which may not be 
practical. On the other hand, we can use a low rovibrational state as \d) state which has less decay channels. Then, 
infrared (IR) spontaneous emission will be the accumulation mechanism instead of optical spontaneous emission and 
the radiation entropy would still increase but with less populated \acc) states and therefore less internal entropy. 
However, the long IR decay lifetime leads to longer cooling time. 

The model presented above does not consider further IR decays from the \acc) states to lower rovibrational states. 
Within about 0.1s, the molecules essentially end up in the ground vibrational state by multiple IR decays. Thus, 
the vibrational heating occuring in the model is a transient effect. We have shown the model up to a stage where 
translational cooling has been achieved to a large extend, which is the aim of our cooling scheme. To summarize, the 
resulting internal heating is not a practical concern in this work since our aim is translational cooling. Translationally 
cooled molecules (even if internally hot) are invaluable assets in many fundamental and applied sciences. 

IX. CONCLUSIONS 

We have developed a laser cooling concept for translational motion of neutral molecules based on the repetition 
of three steps: velocity selection, deceleration and accumulation by a single spontaneous emission per molecule. 
Each molecule need not undergo repeated spontaneous emissions to be translationally cooled. Using realistic cooling 
parameters, we have estimated the cooling time, which turns out to be practically feasible to cool down thermal 
molecules of IK to around recoil width. We have modelled the cooling process and studied the changes in the 
entropies. Simulations of the model have shown the reduction of the momentum width and the center of mass entropy 
by the cooling process. The many internal states of molecules can serve as a useful entropy sink, in addition to the 
radiation entropy for translational cooling. Analysis has shown that our scheme to translationally cool molecules 
by lasers, can be realized in pinciple and in practice. This cooling concept should open up a new possibility for 
translational laser cooling of molecules. 
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FIG. 1. A rotational pumping scheme for initial state preparation 
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FIG. 2. Schematic diagrams for each cooling cycle composed of 3 sequential steps: a)narrow velocity width selection, 
b) deceleration of selected molecules and c) accumulation of decelerated and narrow momentum width molecules by a single 
spontaneous emission to one of the many accumulation states. 
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FIG. 3. Momentum distributions of the internal populations in a typical cooling cycle after the steps of: a)velocity selection, 
b) deceleration and c)accumulation. 
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FIG. 4. Evolution of the a)center of mass entropy S cm , b)internal entropy Si, and c)total entropy Stot as well as S cm ± Si 
throughout the first 2 cooling cycles with 1 accumulation state(line with dots) and 10 accumulation states (solid line). The 
abscissa of each cooling cycle is equally divided for the 3 cooling steps, each with different timescale. 
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FIG. 5. Results after the whole cooling process with 10 accumulation states for a,)S cm and Si, b)Stot and S cm ± Si. 
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FIG. 6. Total momentum probability distribution f(P) before and after the cooling process. 
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